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Abstract 

The MSSM has flat directions in its scalar potential, along which it is natural 
for Bose condensates of squarks to form in the early Universe. A baryon asym- 
metry can be induced in these condensates via Affleck-Dine baryogenesis. The 
condensates are unstable with respect to fragmentation to "B-balls", solitons 
made of squarks and carrying baryon number, which, if they survive thermal- 
ization, fill the Universe down to low temperatures, much lower than that of 
the electroweak phase transition, with interesting cosmological consequences. In 
particular, their decay implies a similar number density of baryons and dark 
matter neutralinos, in accordance with observations. Evasion of thermalization 
and the ability to account for the observed baryon asymmetry requires a very 
low reheating temperature, which, it is argued, is a natural feature of presently 
favoured D-term inflation models. 
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The possibility of electroweak baryogenesis [1]] in the Minimal SUSY Standard 
Model (MSSM) || has become increasingly tightly constrained by experiment. In 
particular, the requirement that the electroweak phase transition is strong enough to 
prevent subsequent wash-out of the asymmetry imposes an upper bound on the Higgs 
mass of 105-107 GeV @ (the present 95 % CL LEP lower bound is 77.5GeV ||, with 
almost all of the remaining range to be tested by LEP200 0), and in addition requires 
that the right-hand stop is light ||, which in turn requires a negative SUSY breaking 
mass squared term for the r.h. stop (but not for the other squarks and sleptons); a 
stiff challenge for supergravity (SUGRA) models. However, there is another perfectly 
natural way to generate the baryon asymmetry in the context of the MSSM, namely 
Affleck-Dine baryogenesis 0. In this the B asymmetry is induced in a squark Bose 
condensate which subsequently decays to (or is thermalized to) quarks. 

This mechanism has recently been given a new spin |7|, |9|] with the realization 
that in almost all realistic cases the squark condensate is unstable with respect to 
fragmentation to B-balls (B carrying Q-balls flit], |IIJ) made of squarks, which fill the 
Universe for a significant part of its history (until well after the electroweak phase 



transition in the case of gravity-mediated SUSY breaking || [L2| and right up to 
and including NOW, with interesting experimental and astrophysical consequences, 
in the case of gauge- mediated SUSY breaking [g, These B-balls can have very 

significant consequences for cosmology; in particular, in the case of gravity-mediated 
SUSY breaking on which we focus here, they can explain the remarkable similarity 
of the number density of baryons and dark matter particles for the case where the 
dark matter particles have masses 0(my/) Hl2| , |i~4| , such as is typically the case for 
neutralinos. To be precise, with currently favoured values compatible with distant 
supernovae observations [13], Q — 1, Q m = 0.4, and Q\ — 0.6, expansion rate h in 
the range 0.6 to 0.8 (as determined by Hubble Space Telescope |16], [17| and limits on 
the age of the Universe |18[) and baryon number density 0.0048 < Qsh 2 < 0.013, as 
determined by primordial nucleosynthesis ITS , we find 

HB --1.5-7.3)^. (1) 
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Thus for tudm ~ m w, the baryon and dark matter particle number densities are 
within an order of magnitude of each other; a remarkable result and a possible clue 
as to the origin of baryons and dark matter. Although one might expect that the 
existence of a significant mass density of baryons could be understood by appealing 
to anthropic selection, there is still no reason for such a selection to result in similar 
number densities of baryons and dark matter. We refer to the scenario in which the B 
asymmetry at least partly comes from late-decaying B-balls as B-ball baryogenesis. 

The details of this mechanism are intimately tied to the details of the full su- 
pergravity (SUGRA) model and inflation, requiring a rather global view of SUSY 
cosmology. This connection with SUSY cosmology as a whole, and the possibility of 
observable predictions |Lj|, |2(| (as discussed in more detail in these Proceedings by 
Enqvist), make this an exciting scenario, requiring, we emphasize, nothing more than 
the particles of the MSSM and a model for inflation. 

B-ball baryogenesis is based on the Affleck-Dine mechanism. This in turn depends 
on the existence of renormalizable flat directions in the scalar potential of the MSSM 



21]. The scalar potential of the MSSM is complicated, with squarks and sleptons 



as well as Higgs fields. Renormalizable flat directions correspond to directions in the 
scalar potential with renormalizible F- and D-terms vanishing, leaving only the soft 
SUSY-breaking terms and the Planck-suppressed non-renormalizable terms which are 
expected in a SUGRA effective theory. These directions are characterized by the 
lowest dimension superpotential operators which are consistent with all symmetries 
and which which are non- vanishing when expressed as scalars along the flat directions. 
We concentrate on the case of the MSSM with R-parity, which eliminates dangerous 
renormalizable B- and L-violating terms from the superpotential and also allows for 



stable neutralino cold dark matter [El 22 . In this case the flat directions are of even 



dimension [21 



e.g. 

d=4 i^L-direction 

<u L >=<<&> ; (H u L) 2 ^0, (2) 



2 



d=6 M c <i c (i c -direction (the d c QL and e c LL directions are similar) 

< u c 1 >=< d c 2 >=< d c ' 3 > ; {u c d c d c f ^ , 



(3) 



where the indices denote colour and d c and d c> are orthogonal in flavour space. We will 
be particularly interested in the d=6 u c d c d c direction, as d=4 directions are both in- 



compatible with B-balls which survive thermalization [12] and with the AD mechanism 



in the context of D-term inflation |23| 



The form of the soft SUSY-breaking terms (once Planck-suppressed SUGRA cou- 
plings are taken into account) is different in the early Universe; the energy density 
of the Universe breaks SUSY and results in 0(H) corrections to the SUSY-breaking 
terms |24|, |25 |. As a result, it is natural for the initial value of the AD squark field to 



take a large value during and/or soon after inflation. This is true if the correction to 
the SUSY breaking mass squared term is negative. The scalar potential then has the 
form 
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v = K - ciJ 2 )!*! 2 + + h.cj + -^l- , (4) 

where M = Mpi/\/8ii and c ~ 1. The AD field will begin to oscillate coherently once 
H « m, corresponding to a Bose condensate of squarks, with an initial amplitude 
(f) « (mM 3 ) 1/4 « 10 14 GeV for the d=6 direction and <\> « (ml) 1 ' 2 w 10 10 GeV along 
the d=4 direction. (It is the larger amplitude along the d=6 direction which protects 
the AD field from being thermalized in D-term inflation models, by giving a large mass 
to the particles which couple to it, so suppressing the scattering rate with light thermal 
particles p3| ). In the absence of the A-terms, there is no B violation or CP violation, 
and the real and imaginary parts of the AD scalar oscillate in phase, corresponding 
to zero B asymmetry. A B asymmetry is induced in the squark condensate if a phase 
shift is induced between the real and imaginary parts of the AD field. The baryon 
asymmetry density is given by 

With $ = ((pi + i(p 2 ) I 'y/2, 4>\ = 4> Cos(mt) and (f) 2 = (f> Cos(mt + 5), the B asymmetry 
density is given by = [m(f) 2 Sin5) jl. The required phase shift is provided by 
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the A-term. When the AD field starts oscillating, the A-term is of the same order 
of magnitude as the mass term and non-renormalizable term in the potential, and 
distinguishes between the real and imaginary directions (where the real direction is 
defined by the phase of the A-term). For the case of interest to us here, D-term 
inflation |26[, in which there are no 0(H) corrections to the A-terms ||23|| , the initial 



phase of the AD field is determined by its random initial value during inflation and is 
therefore typically of the order of 1; thus within a few oscillations a phase difference 
of order 1 will be induced in the AD field. Subsequently, as the Universe expands, the 
amplitude of oscillation decreases; the A-term, which is proportional to 6 , rapidly 
decreases relative to the mass term, effectively switching off B violation and leaving a 
fixed B asymmetry which we observe today. The final B-asymmetry depends crucially 
on the reheating temperature after inflation; for the d=6 u c d c d c direction the baryon 
to entropy ratio is |27| 

Vb « 3 x 10- 11 ( r ^) Sin6 . (6) 

Comparing with the observed B asymmetry, t]b obs = (3 — 8) x 1CT 11 , we see that 
the reheating temperature for d=6 AD baryogenesis must of the order of lGeV in 
D-term inflation models (or in any model in which the magnitude of the CP violating 
phase is of the order of 1; a natural value); for higher reheating temperatures there is 
overproduction of B. 

Thus d=6 AD baryogenesis with 0(1) CP violating phase (as in the D-term in- 
flation scenario) imposes the constraint on inflation models that the reheating tem- 
perature must be very low relative to the energy scale of inflation. At first sight this 
seems perhaps unfavourable; however, reheating in SUSY inflation models must al- 
ready occur at a low temperature relative to the inflation energy scale, T R < IQpGeV , 



in order to avoid thermally overproducing gravitinos |T9|, |28|]. We will see that there 
is good reason to expect a very low reheating temperature in D-term inflation models, 
resulting in an interesting self-consistency of the D-term inflation/AD baryogenesis 



scenario [^J 



In the original AD mechanism it was assumed that the AD field would simply 
oscillate around the minimum of its potential with decreasing amplitude until it either 



decayed or the AD field strength was small enough to allow it to thermalize. However, 
it was recently realized that the AD condensate is generally unstable with respect to 
fragmentation to B-balls || ||. In general, Q-ball solutions (charged non-topological 
solitons in a scalar theory with an approximate U(l) g i b a i, which in our case corresponds 
to baryon number) exist if £7(101 )/|0| 2 has a global minimum at <fi ^ |K|. This is 
likely to occur rather generally in SUSY models, as a result of the soft SUSY breaking 



terms |7j. The resulting solution has the form |T0|| 



$(r,t) = ^=0(r) ; 0(r) « 0(O)e"^ 2 , (7) 
where u is the effective mass of the scalars inside the Q-ball and 0(r) here corresponds 



to the case of thick- walled Q-balls |12 |. For the particular case of flat directions in 
models with gravity-mediated SUSY breaking, the mass squared term rather generally 
becomes smaller with increasing <p once radiative corrections are taken into account, 
as a result of gauge loops [] The RG equations generally have the form |2| 

dim ' 

li-^ = aim] - (5 a M 2 a , (8) 

where rrii represents the scalar soft SUSY breaking terms and M a the gaugino masses. 
These are typically dominated by the gaugino masses, causing to decrease with in- 
creasing scale. Thus the AD scalar mass term decreases with increasing <p. As a result, 
the condition for the existence of Q-ball solutions is satisfied. The homogeneous con- 
densate will naturally fragment to Q-balls, as the Q-ball is the lowest energy solution 
for a given charge. This can be also seen by noting that a homogeneous scalar field 
oscillating in a potential that is flatter than 2 will behave as matter with a negative 
pressure PD| . At 1-loop the correction to the mass term is 

V{\<$>\) = m 2 |$| 2 + ... -> m 2 U + Klog ^rJJ l $ T + ••• > ( 9 ) 

where typically \K\ w 0.01 — 0.1 and K < 0. For small \K\, V ~ <p 2+K , which results 

in an equation of state for the homogeneous field 

K , > 

P = yP • (10) 

x Lct us note some interesing work in connection with flat directions, proton stability, Q-balls and 
string models in this context in p9] . 



This negative pressure merely reflects the attractive interaction between the scalars 
due to the log term in the potential. As a result, any spatial perturbations will grow 
exponentially with time until they go non-linear and the condensate will fragment into 
charged condensate lumps which eventually condense into B-balls. The exponential 
growth of the linear perturbations is given by [|12j 



where 5<p is the seed perturbation, expected to come from quantum fluctuations during 
inflation ||. This is valid for perturbations with |k 2 /a 2 | < \2Km 2 \; shorter wavelength 
perturbations have a positive pressure due to their gradient energy which overcomes the 
negative pressure, such that the perturbations merely oscillate. The first perturbation 
to go non-linear has a diameter l n0 n-Hnear ~ 7r/(2|i ; C|m 2 ) 1 / 2 . The baryon number 
density at a given value of H is given by % = r]BH 2 Mp l /27rTji. From these we find 
that the baryon number of the condensate lumps (and so of the B-balls) is [|Bj 



Z^10 2 VbW /2 (^) , (12) 

with the corresponding radius and field inside the thick-walled B-ball given by R m 
(iK^m)- 1 and 0(0) « 10 14 (fi/10 26 ) 172 GeV. Here f B is fraction of the B asymmetry 
which ends up inside the B-balls; this requires an analysis of the non-linear evolution of 



the condensate [3lJ but we expect B-ball formation to be quite efficient; /# ~ 0.1 — 1. 

Due to their large charges and field strengths, the B-balls can survive thermaliza- 
tion and decay at low temperatures. It may be shown that they survive thermaliza- 
tion if Tr < 10 3-5 GeV [12]. This is easily satified for the d=6 D-term inflation case. 
(However, it rules out survival of B-balls for the d=4 directions, since in this case 
T R > 10 7 GeV is required to account for the observed B [T2"||). Charge escapes from 



B-balls to light B-carrying fermions and scalars a rate proportional to their area, A. 
The decay rate of the B-balls is given by |TT[ [12 



ldB^ cu 3 f s A 1 

B ~dt ~ "192^5 K B ' { ' 
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where & ~ (1 + 3K/2)m « m, with m ~ lOOGeK being the mass of the AD scalar 
[EEH- {fs is a factor accounting for the enhancement of the decay rate if decay to scalar 
pairs is possible; we estimate f s ~ 10 3 Jl2[)- Therefore larger charge B-balls decay 
later. The B-balls decay once |^^| ^ H; the resulting decay temperature is 

r ^(° 01 - 1 )(r^v) I/2GeV - < 14 > 

Thus, in general, B-balls which can evade thermalization will decay after the elec- 
troweak phase transition. This implies that the B-balls will have non-trivial conse- 
quences for cosmology: 

• The baryon asymmetry inside the B-balls (where thermal particles cannot penetrate) 
will be protected from washout due the combined effects of sphaleron B+L violation 
and additional L violating interactions. Such interactions commonly arise in exten- 
sions of the MSSM; for example, in see-saw models of Majorana neutrino masses. 

• ub ~ fiDM- If the reheating temperature is of the order of lGeV, then the B-balls 
will decay after the dark matter neutralinos have frozen out of chemical equilibrium. 
Neutralinos freeze out at Tf m x /20 > 1 GeV, where m x is the neutralino mass p^]. 
The B-ball is made of squarks, with one unit of R-parity per 1/3 B. Thus when the 
B-ball decays we will obtain 3 neutralino LSPs per unit B. This naturally results in 
Ub ~ n-DMi in accordance with observations! Note that this produces non-thermal 
neutralino dark matter, with a quite different dependence on the MSSM parameter 
space than the conventional thermal relic dark matter, a feature that can be tested 
experimentally once SUSY is discovered in accelerators. The actual number density 
ratio predicted by B-ball decay (assuming no subsequent annihilations |12| ) is crucially 
dependent on f B 

- = tt- < 15) 

n x 3f B 

Thus fs has to be less than 1 in order to be compatible with observations (Eq.(l)) and 
the present experimental lower bound on the LSP neutralino mass, m x > 25GeV r [5^] . 
However, fs in the range 0.1 to 1 can accomodate a wide range of neutralino masses. 
In addition, if we can calculate fs, we can constrain the LSP mass. For example, for 



the Vt m = 0.4 case we find 

771 

0.046/s 1 < < 0.22/^ , (16) 
raw 

where the range is due to the uncertainties in the baryon asymmetry estimate from 
nucleosynthesis and in the present expansion rate. Since Q m , Qb and h will all be 
fixed to 1 % accuracy by the MAP and PLANCK satellite missions [[33], |IT|, we will 



eventually be able to predict m x if we can calculate providing a " smoking-gun" 
test of this scenario |l4j . It remains to be seen whether fg can be accurately computed 



from the non-linear evolution of the condensate . 

• Another observable consequence of late-decaying B-balls, in the context of D-term 
inflation models, is that they transfer the isocurvature fluctuations in the baryons to 
the dark matter neutralinos, resulting in a significant enhancement of the isocurvature 



fluctuations, making them observable by MAP and PLANCK ||20|| . This is discussed 
in detail in these Proceedings by Enqvist. 

So far we have shown that the d=6 Affleck-Dine mechanism with 0(1) CP violating 
phase requires a reheating temperature of the order of lGeV, and we have noted that 
0(1) phases are to be expected in D-term inflation models. We now address the 
question of whether such low reheating temperatures are a natural feature of these 
models. 

SUSY inflation is characterized by the nature of the energy density driving inflation. 
The SUSY potential has F- and D-term contributions, V susy = \F\ 2 + \D\ 2 , with D- 
term inflation driven by \D\ 2 . The reason D-term inflation is favoured over F-term is 



that there are no 0(H) corrections to the inflaton mass in this case [26|. These terms 
would otherwise ruin the flatness of the inflaton potential, preventing slow-rolling 
and/or producing too large deviations from a scale-invariant perturbation spectrum. 

The minimal model for D-term inflation has a singlet S and two fields with opposite 
U(1)fi charges ip±, where U(1)fi is the Fayet-Illiopoulos (FI) gauge group with a FI 
term £ 2 . The superpotential is 

W = AS^+V- (17) 



s 



and the scalar potential, including the U(1)fi D-term, is then 

V = |A| 2 (|^_| 2 + |^+| 2 + |S^| 2 ) + ^(|^ + | 2 - |^| 2 + e? , (18) 

where g is the U(1)fi coupling. £ is fixed by cosmic microwave background measure- 
ments to be 6.6 x 10 15 GeV [3*3] . The inflation model which results from this potential 



has the form of a hybrid inflation model ||35|| ; for S > S cr u = g£/\ the minimum of 



the potential corresponds to ip± = 0, resulting in a very flat, slow-roling potential 
with V = g 2 £ 4 /2 + k log\S\, where k comes from radiative corrections. For sufficient 
inflation, S > S55 ~ 0.9M is necessary. However, this leads to a generic problem for 
D-term inflation models . In general, since we are considering SUSY models to be 
the low energy effective theory of a SUGRA theory, we expect all non-renormalizable 
superpotential terms consistent with the symmetries and suppressed by powers of 
M to appear, AW ~ S m /M m ~ 3 . However, the contribution to the scalar potential 
of such terms will cause unacceptable deviation from scale-invariant fluctuations at 
length scales corresponding to the observable Universe, requiring the elimination of 
all such terms for m < 9 fl23 |. Thus we need to introduce a symmetry to eliminate 



such terms; we focus on the case of an R-symmetry, which is particularly effective 
(eliminating all purely S superpotential terms if R(S) is negative, for example), al- 
though we expect the same argument to apply for more general discrete symmetries. 
The connection between D-term inflation and low reheating temperatures is then that 
the symmetry which keeps the inflaton potential sufficiently flat also tends to suppress 
the coupling of the inflaton to light fields, so suppressing the decay rate of the infla- 
ton and resulting in very low reheating temperatures P71 []. This happens if there is a 
"mismatch" between the transformation properties of the inflaton sector S, ip± and the 
MSSM sector fields. As a simple example, consider the case where the inflaton sector 
fields have half- integral charges and the MSSM sector fields have 1/3-integral charges. 

2 In D-term inflation models there are two stages of reheating; that from decay of the tp- oscilla- 
tions, with typically Tjj w 10 15 GeV, and that from decay of the S oscillations p3| . The first stage 
of reheating is important for thermalizing the d=4 AD field, but otherwise the important reheating 
temperature (for the gravitino bound and for B-balls) is that associated with the S field, whose 
oscillations dominate the energy density throughout p7|] . 
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With R(S) = -n, R(ip + ) = R(if>-) = (n + 2)/2 (n e Z) for the inflaton sector fields 
and R(Q) = R(L) = 5/3, R(u c ) = R(d c ) = R(e c ) = 1/3 and R(H U ) = R(H d ) = 
for the MSSM fields (the MSSM charges have been chosen to allow the (u c d c d c ) 2 op- 
erator necessary for AD baryogenesis) we find that, for n=l, the lowest dimensional 
superpotential term allowing the inflaton to couple to the MSSM sector is of the form 
Wint — nS(LH u )(u c u c d c e c ) / M A ~ nScff / M r ~ 2 with r=6, where represents the light 
MSSM fields |[27|| . For n=2 a dimension r = 4 operator is possible, but for larger n the 
dimension is r > 6. The corresponding decay rate is 

{M s \ 2 ^ 



s 



\ M 



K 2 (3 r M s , (19) 



where (3 r represents the phase factor for decay to r light particles. The corresponding 
reheating temperature is then 

T| w 4 x 10 15 (2.8 x l0-3) r - 5/2 K\ r - 3/2 (3l /2 GeV . (20) 
Thus for r=6 we find p7 | 



/ A \ 9/2 

Tl « 150 -j K /? 6 1/2 GeV , (21) 
where (3q ~ 10~ 6 . Therefore low reheating temperatures, as low as 1 GeV or less, are 



a natural feature of D-term inflation models |27] . 

Of course, it is also necessary to check if reheating mechanisms other than sin- 
gle particle decay could be more efficient; in particular parametric resonance via 
Planck-suppressed operators. We find that parametric resonance is completely negli- 
gible if A/g 1 / 2 ^ 0.03 PT| . For values of A and g not too small compared with 1 this 
condition will be easily satisfied. 

In conclusion, we believe that Affleck-Dine baryogenesis with late decaying B-balls 
(B-ball baryogenesis) is a natural candidate for baryogenesis in the MSSM which has 
the great advantage that it can explain why ~ t^dm- This model fits in very 
well with the currently favoured D-term inflation scenario, with a natural consistency 
between the low reheating temperature required for d=6 AD baryogenesis and the 
suppression of inflaton couplings required for a flat inflaton potential. It is especially 

10 



interesting that this model can make definite predictions; % ~ udm, non-thermal 
neutralino dark matter, eventually (depending on the calculation of /# and results 
from MAP and PLANCK) the neutralino mass, and, as discussed in these Proceedings 
by Enqvist, isocurvature microwave fluctuations which should be observable by MAP 
and PLANCK. The correlation between these various predictions should serve as a 
clear test of the B-ball baryogenesis/D-term inflation scenario. 
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